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ABSTRACT 


Spiraling in ion-molecule collisions with polarizable and/or polar targets is inves- 
tigated numerically using computer plotting techniques. Collision orbits and interaction 
behavior are studied for a range of collision parameters on the IBM 360/67 and 7094 
computers with the CDC DD280 plotter. Multiple reflections off a hard-core barrier at 

small ion-molecule separations occur only in collisions with permanent dipoles. 

3 

Multiple- reflection times for 2 to 2000 reflections are 10 to 10 times longer than times 
for single reflections. Long multiple-reflection collision times, however, do not mean 
large capture cross sections; the capture cross-section values are much less than those 
predicted by the adiabatic theory. 



A COMPUTER PLOTTING DESCRIPTION OF ION-MOLECULE COLLISIONS 
WITH LONG-LIVED CAPTURE COMPLEXES 
by John V. Dugan, Jr., and James H. Rice 
Lewis Research Center 

SUMMARY 

The degree to which spiraling occurs in ion- molecule collisions with polarizable 
and/or polar targets is investigated numerically by means of computer plotting tech- 
niques. The collision orbits and ion-molecule interaction behavior are studied for a 
range of collision parameters by using the IBM 360/67 and 7094 computers with the 
CDC DD280 plotter. The interaction potential consists of a hard-sphere term plus the 
attractive polarization and ion- dipole terms. Multiple reflections off the hard-core 
barrier at small ion-molecule separations (i. e. , r « the reflection distance r c ) occur 
only in collisions with permanent dipoles. Simple specular reflection occurs in pure 
polarizability (Langevin) collisions. No spiraling, as conventionally defined, is ob- 
served in the approximately 700 collisions studied. Multiple- reflection times t are 

3 ^ 

10 to 10 times longer than the corresponding values for single reflections. As many 

as 2000 reflections are observed at r c values of 1 to 3 A. Long multiple- reflection 

collision times, however, do not mean large capture cross sections o c ; the capture 

cross-section values are much less than predicted by the adiabatic theory. 


INTRODUCTION 

Collisions between ions and polarizable molecules have been considered as spiraling 
encounters since the early paper by Langevin (ref. 1). This view was formulated (ref. 2) 
from the simple orbit equation without any discussion of the conditions necessary for 
such spiraling behavior. Results of numerical calculations of capture cross sections 
for ion - permanent dipole collisions were compared (ref. 3) with cross sections for 
Langevin (dipole moment, p = 0) collisions involving hydrochloric acid (HC1) targets. 

The nature of these collision orbits was also examined. These results showed no sig- 
nificant spiraling for the Langevin test cases in the absence of a repulsive potential. 


A recent proposal (ref. 4) of accidental resonance in ion-molecule collisions sug- 
gests spiraling behavior to explain very long collision times for charge exchange. This 
report provides results of a numerical investigation of ion-molecule collisions by means 
of plotting techniques. Collisions involve both polar and nonpolar molecules. The ex- 
tent of spiraling is investigated, as well as conditions for reflection of the particles off 
a repulsive barrier. The nature of the collision complex formed in multiple- reflection 
captures is treated in some detail. A cursory study of the effects of target geometry on 
ion-molecule complex formation is also included. All symbols are defined in appen- 
dix A. The SI system of units is used throughout except for length where angstroms 
(1 A = 10" 10 m) and energy where electron volts (1 eV = 1.602x10”^ J) are employed 
because of the atomic nature of the problem. 


COLLISION MODELS 
Langevln Collisions 

The interaction between an ion and a polarizable molecule (or atom) is described by 
the ion-induced dipole potential 

V(r) “5? — (1) 

2(47re o )r 4 

where a is the electronic polarizability, e is the electronic charge, and r is the ion- 
molecule separation. In the following discussion this potential energy term is written in 
atomic units so the mks conversion factor of (47re o ) _ ^ will be omitted. This attractive 
potential leads to an orbit equation (ref. 2) which relates the scattering angle /3 to the 
separation r by 



(2) 


where p = b/r and p Q is the lower of the two positive roots of the radical in equa- 
tion (2), if the roots exist. For b < b T (where b T is the Langevin critical impact pa- 

2 1/4 " ■ L * 

rameter equal to (2 ae /e) ' ), there is no real root and the orbit is of the spiraling 
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type. Orbits for b > never approach closer than r R = and e is the initial 

relative translational energy of the ion and molecule. Thus the probability of capture 
C R drops off abruptly to zero (as a step function) at b = b R . Gioumousis and Stevenson 
(ref. 5) calculated the cross sections for various ion-molecule reactions by assuming 
that the critical reaction radius is equal to r^. They obtained satisfactory agreement 
with certain experiments. Hamill and his co-workers (refs. 6 and 7) studied the ques- 
tion of reaction probabilities as functions of relative velocity experimentally and ob- 
tained results consistent with the calculations of references 3, 8, and 9. 


Permanent Dipole Collisions 

The ion-dipole contribution to the interaction potential is 

V(r, y) = - cos . y 
47re o r 2 


(3) 


where (i is the dipole moment and y is the angle between the separation vector and the 
negative end of the dipole. The sign of the permanent dipole potential naturally changes 
as the polar molecule rotates. The collision orbits for polar targets are not amenable 
to the simple spiraling description of equation (2). The probability of capture C R for 
such collisions has been calculated as a function of impact parameter at fixed values of 
relative velocity (refs. 3, 8, and 9). This C R value is a complicated function of im- 
pact parameter compared to the step-function dropoff in capture probability for Langevin 
collisions. 

It is clear that the presence of a permanent dipole is necessary for multiple reflec- 
tions to occur. In the absence of a dipole, the behavior of the ion-molecule pair is best 
understood in terms of the effective potential V efP which can be written 


v eff< r > = 


-er-i 


ae 


(4) 


The maximum in this potential is at 


r . 
mm 



(5) 
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Capture (defined as approach to r < r c in this report) occurs only for b < b^; however, 
for these trajectories the particles simply reflect and pass back over the barrier at 
r = r*. The condition for capture for b < b^ is, in fact, that the relative energy e 
T 2 / 41^/2 

be greater than 2V eff (r*)ofe /b J .No spiraling is observed in collisions without di- 
poles (ref. 10) (i.e., pure polarizability collisions). For such collisions (p = 0), it is 
concluded in reference 10 that long collision times in charge exchange collisions cannot 
be explained in terms of spiraling behavior. 

The presence of a rotating dipole considerably complicates the effective potential 
description of the ion- molecule collision. Multiple maxima can be introduced in the 
effective potential; the location of these maxima changes as the dipole rotates. The 
height of the potential barrier can shift as well. This rapid time variation of the poten- 
tial permits trapping of the ion after reflection off the step potential (i. e. , multiple re- 
flection can occur). A barrier in the effective potential can be introduced for any 
r =- r . 


Collision Times 

The question of the duration of collisions has arisen in theoretical predictions of 
Bohme, Hasted, and Ong for cross sections in charge exchange processes (ref. 4). 

These reactions include polar and nonpolar targets. It is assumed in reference 4 that 
the adiabatic theory is a satisfactory description of these collisions. "Adiabatic” in 
this sense means that the electronic configurations of ion and molecule can adjust as the 
collision proceeds. They then employ Massey’s adiabatic description of the cross sec- 
tion, which says that it reaches a maximum when the collision time t q equals the char- 
acteristic transition time t^. The authors of reference 4 suggest that ion- molecule 
charge exchange proceeds as an "accidentally resonant” process because spiraling im- 
plies long collision times. In the adiabatic approximation, long characteristic transition 
times imply small transition energies AE^ » 0 (i. e. , resonant collisions). (Since the 
resonance is not due. to symmetry effects, the collisions are termed "accidentally reso- 
nant. ”) 

Five of the reactions considered in reference 4 involve the polar target molecules 
carbon monoxide (CO) and nitric oxide (NO). Both of these targets have relatively small 
dipole moments, 0. 1 and 0. 16 Debye units. The present study considers cases of low 
(CO), intermediate (HC1), and high (methyl cyanide (CH^CN)) dipole moments. The case 
of CO can be compared with the predictions of reference 4 for charge exchange cross 
sections. The N02 + - HC1 system was chosen not only as representative of intermediate 
H value but because the reduced mass for this pair is identical to the value for the 
CH 3 CN + + CHgCN case. 
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Repulsive Potential 


The numerical trajectory calculations of references 3, 8, and 9 were performed 
with the induced dipole potential term, permanent dipole term, or both terms. If in the 
course of the trajectory, pair separations of 1 to 2 A are achieved, the effects of a re- 
pulsive potential become significant. The electronic clouds of ion and molecule overlap 
at these r values. In these previous studies (refs. 3, 8, and 9), the trajectories were 
stopped where r = r c (ion-molecule separations of 1 to 2 A where repulsive effects oc- 
cur). In the present study, a hard-sphere potential barrier is simulated by changing the 
sign of the radial velocity r at r = r . Spiraling is defined as changes in the polar 
angle 0 by v and/or the azimuthal angle cp by 2tt. In this report we wish to clearly 
distinguish spiraling (which does not involve multiple reflection) from a trajectory where 
multiple reflections result in a circuit of the ion around the molecule. 


NUMERICAL TRAJECTORY CALCULATION 

Details of the numerical calculation have been described in references 3, 8, and 9. 
The Lagrangian equations of motion for an ion interacting with a linear polar molecule 
are given in appendix B. The effects of geometry on capture cross sections for ion- 
dipole collisions have been discussed in references 8 and 9; the equations of motion are 
very similar for models of symmetric top molecules (ref. 9). The coordinate systems 
used in the computer study are shown in figures 1(a) and (b). Molecular constants and 
collision parameters are given in table I. 


TABLE I. - COLLISION PARAMETERS 1 ’ b 


Ion- molecule pair 

Electronic 

polariz- 

ability, 

a, 

(A) 3 

Dipole 

moment, 

4. 

Debye unit 

Reduced mass, 
m, 
kg 

Moment of inertia, 

I, 

kg-m 2 

Ar + + CO 

1.95 

0. 10 

2.74X10' 26 

1.44xl0' 46 

N0 2 + + HC1 

2.60 

1.08 

3.41 

.268 

CH 3 CN + + CHgCN 

3.80 

3.92 

3.41 

9. 12 


a All collisions were studied for an initial ion-molecule relative velocity 
v «5xl0 4 cm/sec; the reflection barrier was located at r = 1, 1.5, 

2° or 3 A. 

^Target rotators were chosen from a distribution at a rotational tempera- 
ture Tpj = 300 K. 
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Trajectories were started at r Q = 50 or 25 A and continued after interaction until 
r r Q . If the ion-molecule, pair reaches a turning point for r > r c (the "capture" dis- 
tance), the collision is recorded as a "repulsion. " Orbits for which r < r c at any 
time are denoted "capture" collisions. The step-size changing routine (variable order 
Runge-Kutta technique) used to integrate the differential equations is described briefly 
in appendix C. 

The collision time t q is also of interest for these interactions since according to 
the adiabatic theorem (refs. 4 and 11) large r c implies large reaction cross sections. 
The choice of interaction distance r* (with t q defined as the period for which r < r*) 
is somewhat arbitrary because of the presence of the dipole. It will be shown, however, 
that for multiple- reflection collisions the t c value is relatively insensitive to the choice 
of r*. A lower limit to r* (r ^- n of eq. (5)) is the maximum in the one- dimensional ef- 
fective potential for pure polarizability collisions given by equation (4). This potential 
can be written 


2 T >2 
-ae L 

= + 

eff a c 


2r 


2mr 


where L* is the angular momentum (L' = mvb) and m and v are the reduced mass 
and relative velocity. An alternative choice of r* is that distance at which transla- 
tional velocity is first significantly perturbed (~10 percent). This distance varies be- 
tween 4 to 20 A for the three collision systems studied. Its value does not differ signifi- 
cantly from the hindering distance r h (refs. 3 and 8) for representative rotators where 
r^ is defined by 


f “ e r < 6 > 

r h 

As long as r c is large (i.e. , r c » t q , the average collision time for a single reflec- 
tion), the r c value remains insensitive to the choice of r*. Capture collisions with 
t c » t q are not strictly specular reflections since their orbits are not so simple as in 
the pure polarizability (Langevin) case. These capture collisions are designated 
"single- reflection collisions" in the following sections. 


COMPUTER PLOTTING APPROACH 

The collision trajectories were solved on the IBM 360/67 and 7094 computers and 
displayed by using the CDC DD280 plotter. The plots of translational motion are done 
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in a coordinate system where the polar molecule is fixed at the origin. The spherical 
coordinate angles 6 and cp for translational motion were plotted to determine if spiral- 
ing behavior occurred. The relative velocity v, rotational energy E-^, rotator projec- 
tions, orientation angle y, and ion projections in Cartesian planes could all be plotted 
so that each feature of the collision could be scanned. Preliminary work on nonreflected 
trajectories stopped at r = r c was reported in reference 10. The procedure for plotting 
and a discussion of running times for multiple- reflection collisions are included in appen- 
dix C. 


RESULTS AND DISCUSSION 

It is of interest to study the nature of multiple reflections as well as the shape of the 
ion orbits in the field of the dipole. The data for collision times and number of reflec- 
tions in multiple-reflection captures are included in table EL These data are representa- 
tive of all the collisions which were studied for the collision parameters given in table I. 
These results will be discussed in detail after the nature of orbiting and energy exchange 
have been discussed for all three polar targets. 


Collisions of Ar + With CO 

The variations of polar 6 and azimuthal cp angles (for translational motion) with 
ion- molecule separation for an Ar + + CO capture collision are shown in figure 2. The 

4 

relative velocity is 5x10 centimeters per second, corresponding to a relative transla- 
tional energy e « 0.03 electron volt (near thermal). The energy of the rotator in this 
case was «0. 5 kTp^ where Tp^ = 300 K (relatively cold). The initial conditions for all 
variables and time derivatives are as follows: 

Initial conditions: 


r = 25 A 
6 = 1 . 571 rad 
<p = 4. 140 rad 
7] = 5. 870 rad 
£ = 2.511 rad 


b = 6 A 

v = 5.5X10 -2 A/r n 

E-p = 4.99X10" 6 m AVr 2 
R n / n 

m = 1.44X10" 3 m n 
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TABLE II. - RESULTS FOR MULTIPLE-REFLECTION COLLISIONS 


Ion-molecule pair 

Impact 

Location 

j Number ! Langevin 

Range of > 

Fraction 

Maximum 

Average single- 

Most 

Average collision 

Maximum 


param- 

of 

1 of | 

potential 

numerical | 

of cases 

number 

reflection time, 

probable 

time for multiple 

collision 


eter, 

reflecting 

cases, 1 

maximum, 

interaction j 

resulting 

of 

V 

number of 

reflection, 

1 

time, 


b, 

barrier, 

n 

r*, 

radius, j 

in 

reflections, 

sec 

multiple 

'R. 

T max’ 



r c> 

i 

A 

defining j 

multiple 

^max 


reflections, 

sec 

sec 



A 



collision 

reflection. 



N„ 




Ar - CO 



10 

8.5 

~ 


40 

8.5 J 

8 t 


8 

5.7 

~ 


56 

5.7 

~ 


6 

4.3 

~ 


3 

15 

1.5X10 

6 

3 

5 

2.5 

2 


0 

1115 

1.5x10 

20 

0 

40 

2 

4 

0 

192 

2 

7 

4 

530 

1.7 

3 

7 

715 

4 

6 


to 7x10 


15 to 20 
15 to 20 
15 to 20 
15 to 20 
15 to 20 
15 to 20 
15 to 20 


Case limited to 25 reflections to minimize computer time. 
















180 




(b) Azimuthal angle, <p. 

Figure 2. - Variation of coordinate angles for argon ion (Ar + ) relative to 
carbon monoxide (CO) molecule during single-reflection capture 
collision. 
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Time derivatives: 


r = -5. 34xl0" 2 A/r n \ = 4. 25xlO' 2 rad/T n 

6 = 4.65X10" 4 rad/r n r? = 1.63X10" 1 rad/T n 

<p = 2. 5X10 " 4 rad/r n 

14 -20 -23 

where the scaling values are r = 10" second and m = 10" gram (10" kg). The 

‘I 11 Q 

capture trajectory of figure 2 was solved for an r value of 1 A at an impact parameter 

o c —12 

of 6 A. it is a single reflection where the collision time t « 1. 5x10 second. This 

is representative behavior for such collisions where the dipole interaction is similar to 

the polarizability potential for b < b^. The collision distance r* («5 A) was obtained 

from the plot of translational velocity v against separation distance (i.e. , where e is 

perturbed by 10 percent). The collision time calculated numerically is not very different 

-1 2 

from r Q « 2r*/v Q = 10" second since r c « r*. 

Energy exchange between translational and rotational modes can be conveniently dis- 
cussed in terms of the relative velocity and translational energy. Figure 3(a) shows the 
variation of relative velocity v with ion-molecule ser<&»*ation. There is considerable 




Ion-molecule separation, A 
(b) Molecular rotational energy. 


Figure 3. - Variation of relative velocity and molecular rotational energy 
during Ar + + CO single-reflection capture collision. 
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acceleration of the translational motion. This acceleration generally occurs within 5 A. 
The maximum acceleration is determined by the Me/ r value at r « r . Very little 
rotational energy seems to have been transformed to the translational mode as the par- 
ticles return to r . Correspondingly, although the rotational energy is usually consider- 
ably increased during interaction (as shown in fig. 3(b)), the rotor usually returns to its 
initial energy state. This behavior is consistent with the experimental condition of a 
monoenergetic beam of ions with velocity equal to y 2kT R /m incident upon a Maxwellian 
distribution of rotators at T R = 300 K (refs. 6 and 7). The oscillatory pattern of varia- 
tion in rotational energy is characteristic of these ion-dipole collisions. The frequency 
is inversely proportional to the moment of inertia, except at small r where the varia- 
tion becomes quite random. 



-1 o 1 

X axis 




Z axis 



Ion-molecule separation, A 


Figure 4. - Variation of rotator projections and orientation angle during 
Ar + + CO single-reflection capture collision. 


Figure 4 shows the variation in the rotator projections as the collision proceeds; the 
curves are not completely smooth because every second point was plotted. The CO tar- 
get has a moderately large moment of inertia so the rotator trace shows only a modest 
number of rotations during collision. In several collisions there is evidence of hindered 
rotation of the dipole by the incident ion (i. e. , dipole is preferentially alined with its 
negative end facing the ion) . 

Figure 5 shows the interesting orbiting behavior for an Ar + + CO capture collision 
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where the ion-molecule pair is multiply reflected off a stepwise barrier. This particular 

o 

encounter resulted in eight reflections. For this critical distance (r c = 2 A), as many as 

40 or 50 reflections occurred in some cases for CO targets (see table n); however, it is 

difficult to determine the time history from a 6 against r plot with so many reflections. 

The maximum b values at which a significant fraction of captures occur are ap- 

2 

proximately equal to the b values where the numerical cross section a = nb . The 

C r C C 

a values drop off to quite low values for v = 2x10° centimeters per second. There are 

v 

no multiple -reflection capture collisions observed for CO at this velocity. The ion- 
molecule relative energy £ is so large at these small b values that the dipole term is 
prevented from introducing reflecting maxima in the effective potential. The a c values 

are approximately equal to a T since interaction at the turning points is dominated by 

2 4 -1/2 

the -cte /2r term. The slope of against e is roughly e“ ' , however, which is 

similar to the Langevin prediction for the polarizability potential (ref. 1) . A change in 

slope might be expected if there is a maximum amount of translational energy which can 

be randomized (absorbed) as internal energy in the ion-molecule complex. 


Collisions of N0 2 + With HCI 

Orbiting behavior in a multiple -reflection collision for the system N0 2 + + HCI is 
shown in figure 6. The ion -dipole term dominates the interaction here, in contrast to 
the CO case. The traces do not return to r Q because the plot capacity was exceeded 
for this case. 

The general features of the NC> 2 + + HCI interaction are shown in figure 7. The rela- 
tive velocity and rotational energy first become perturbed at 5 to 10 A because the domi- 
2 

nant -pe/r term is significant at these distances. The variations in both v and Ej^ 

are similar to the case of CO targets. The heavy structure in the rotator projection 

-47 

plots of figure 8 is a result of the small moment of inertia for HCI, 2. 68x10 kilogram - 
2 

meter . The motion of the rotator is clearly as random as for the CO case, although 
some hindered rotations are also observed. 

Collisions of CH 3 CN + With CH 3 CN 

Multiple -reflection behavior for CHgCN targets is shown in figure 9 by means of the 
variation of polar and azimuthal angles with ion-molecule separation. The striking char- 
acteristic of these capture collisions is the large values at which the turning points occur 
for multiple reflections. The maximum observed turning-point value is 22 A. The re- 

o 

flections in figures 9(a) and (b) have a maximum turning point at about 10 A. Figure 10 
shows the case of double minima where the large dipole moment introduces turning points 


14 


Polar angle, deg 




Figure 6. - Variation of coordinate for nitrogen dioxide ion (NOj" 1 ") rela- 
tive to hydrochloric acid (HCI) molecule during multiple-reflection 
capture collision. 
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(3) Relative velocity. 



5 20 15 10 0 5 0 

Ion-molecule separation, A 

(b) Molecular rotational energy. 

Figure 7. - Variations of relative velocity and molecular rotational energy 
during N0 2 + + HCI multiple-reflection capture collision. 
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Z axis Ion-molecule separation, A 

Figure 8. - Variations of rotator projections and ion-molecule orienta- 
tion angle during N0 2 + + HCI multiple-reflection capture collision. 














at 2 to 6 A and at 18 A. 

Again the features of the v and E^ variation are similar to other targets; yet be- 
cause of the large p value, interaction begins at large r. Because of the large moment 
of inertia for CHgCN (five times the CO value), there are usually relatively few traces in 
the rotator projection plots. However, multiple reflections obscure the time history, 
and hindered rotation is not readily apparent (see fig. 11). The random (increasing or 
decreasing) amplitude variation in orientation angle y is characteristic of all polar tar- 
gets. Generally, as illustrated in figure 11, the mean value of y for small separations 
is between 0° and 60° (i. e. , the dipole is in an attractive position with respect to the ion.) 



Zaxis Ion-molecule separation, A 

Figure 11. - Variations of rotator projections and ion-molecule orienta- 
tion a,ng!e during CHjCN + + CH 3 CN multiple-reflection capture colli- 
sion with double minima. 

Numerical Collision Time Results 

Table n includes values of collision times for single- and multiple-reflection colli- 
sions, fraction of collisions resulting in multiple reflection, and mean times and reflec- 
tion numbers for multiple -reflection collisions. Each set of cases consisted of the indi- 
cated number of randomly generated cases at a fixed impact parameter. The most prob- 
able number N M and maximum number N max of reflections are also given; 

The fraction of capture collisions in which multiple reflections occurred for the 162 






Ar + + CO cases studied varied from 0. 2 to 0. 8 at impact parameters from 2 to 6 A. 

o 

This fraction falls off abruptly for 6 to 7 A because the capture cross section is limited 

to b sts b T . These results are generally true for r values of 1 and 2 A. For large 

r c values (3 to 4 A), there are no reflections in the range of interesting b values since 

the interaction is dominated by the polarizability potential at the turning points. It should 

be noted that long collision times r » t , which are observed, do not imply large 

c o 

cross sections a . The capture cross section for CO and NOg targets is approximately 
the Langevin cross section independent of the occurrence of multiple reflections. The 
characteristic radii (turning points) of the ion-molecule collision complex for Ar + + CO 
are 3 to 6 A (ref. 10). 

Single -reflection times for HC1 targets are approximately equal to the CO values at 
the same velocities (see table n). However, HC1 target collisions usually show many 
more reflections than CO targets. This is partly due to the small moment of inertia for 
HC1 (i. e. , high rotational frequency). Also the high dipole moment causes a stronger 
interaction with larger changes in value. Table n allows for a comparison of colli- 
sion times for multiple and single reflections involving CO (ju = 0. 1 Debye unit) and HC1 
(/i = 1.08 Debye unit). It is clear that the fU values are generally larger for HC1 (com- 

o ^ l 

pared to CO) at r fi values of 1 and 2 A. The characteristic radius r t of the HC1 + N0 2 
complex, «5 A, is slightly larger than for Ar + + CO because the higher dipole moment of 
HC1 introduces potential maxima at larger r values. The variation in ion velocity is 
such that the collision time cannot be simply computed from the number of reflections 
with the initial velocity value. 

It is not clear why the multiple -reflection behavior for HC1 is strongly dependent on 
r value. Also interesting but not explainable for the NO 0 + + HC1 collisions is the very 
large number of multiple reflections that occur for some cases; it appears that, once 
trapped, the ion finds it very difficult to escape from the HC1 molecule. 

The t values for multiple -reflection CH„CN collisions are generally large com- 
pared to single -reflection times but not quite as large as the corresponding r c values 
for HC1 targets. The fraction of capture collisions which result in specular reflection is 

_ ' o 

relatively large for CHoCN collisions with r values of 1 and 2 A and remains signifi- 

o ^ ^ 

cant even for r c = 3 A. It is noteworthy that this complex has a characteristic radius 
which is considerably larger than molecular bond distances. It would certainly be of in- 
terest to calculate charge transfer probabilities for these collisions since the heavy- 


particle motion has been calculated. The capture cross section for CHgCN targets 


(ref. 9) is approximately equal to the area described by the orbit radius 
15 A. These capture cross-section values (& c = nb^j are approximately 


b fi of 12 to 
1000 square ang- 


stroms at thermal energy. 
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CONCLUDING REMARKS 


The approximation of a classical ion-dipole plus ion-induced dipole (polarizability) 
potential has been used with a hard -sphere reflecting barrier to study collisions. With 
this approach, it has been demonstrated that spiraling (as distinguished from multiple 
reflections) of collision partners does not occur over the wide range of collision param- 
eters studied. Multiple reflections off a repulsive barrier do occur but only for colli- 
sions with polar targets. The fraction of multiple -reflection collisions depends on the 
hard-sphere radius r £ and the electrical properties of the target atom (i.e. , polariza- 
bility and dipole moment). The CO molecule represents the limit of small /j (0. 10 
Debye unit) when the polarizability term dominates the interactions. The capture cross 
section in this case is approximately Langevin, but the dipole still causes multiple re- 

o 

flection to be probable if r £ ^ 2 A. The CHgCN molecule is the large dipole moment 

limit when the ion -dipole term dominates, leading to large c t values (a »a T ). The 

probability of multiple -reflection collisions is larger for CHgCN than for CO as long as 

r <3A for the former. Collision times for multiple -reflection collisions are longest 
c 3 

for HC1, as much as 10 times single -reflection times, but the fraction of such collisions 
for HC1 targets is small. The characteristic radii of ion-molecule complexes are large 
(3 to 22 A) and generally approximate the critical impact parameters for capture. The 
presence of the dipole is clearly required for multiple reflection with the mathematical 
model used. The general features of the collisions (i.e. , variation in relative velocity 
and rotational energy) are similar for all polar targets. 

Lewis Research Center, 

National Aeronautics and Space Administration, 

Cleveland, Ohio, June 5, 1969, 

120-26-03-04-22. 
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APPENDIX A 


SYMBOLS 

o 

b impact parameter, A 

b c critical impact parameter defining numerical capture cross section, A 

o 

b^ Langevin critical impact parameter, A 

C R capture probability, dimensionless 

-IQ 

E r rotational energy of polar molecule, eV (1 eV = 1.602x10" J) 

AE^. transition energy for ion-molecule charge transfer, J 

1 Q 

e electronic charge, 1.602x10 C 

I o 

F y potential function in translational acceleration expression, A 
k Boltzmann constant, 1. 381x10"^ j/K 

L Lagrangian energy expression, T - V, eV 

o 

L t translational angular momentum, mvb, kg-m /sec 

m reduced mass, kg 

q generalized coordinate 

o 

r ion-molecule separation, A 

o 

r c location of hard-core potential, A 

o 

r^ hindering distance for polar rotation, A 

o 

r-^ distance of closest approach in Langevin collisions, A 

o 

r Q initial ion-molecule separation, A 

o 

r^ turning-point radius for ion-molecule collision complex, A 

• o 

r* numerically computed interaction distance, A 

r min l° ca -tion of maximum in effective polarizability potential, A 
T kinetic energy of ion -molecule system, rotational plus translational, eV 

T r rotational temperature of polar molecule targets, K 
V ion-molecule interaction potential energy, eV 

V e ff effective one -dimensional potential energy of colliding system, eV 
V 0 hard-sphere potential barrier of infinite height 
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V 


relative translational velocity, m/sec 
v Q initial relative translational velocity, m/sec 

03 

a electronic polarizability, (A) 

|6 scattering angle for Langevin collision, rad 

y orientation angle between dipole and ion-molecule radius vector, rad 
e relative translational kinetic energy, eV 

-12 -1 

dielectric permittivity of vacuum, 8.537x10 farad meter 
77 azimuthal angle for rotational motion, rad 
9 polar angle for translational motion, rad 

30 

p. permanent electric dipole moment, C-m (1 Debye unit = 3. 33x10“ C-m) 

| polar angle for rotational motion, rad 

p reduced impact parameter (for Langevin collisions p = b/r), dimensionless 

0 2 

cr cross section, (A) 

o 9 

ct c numerical capture cross section, (A) 

0 2 

a L Langevin capture cross section, (A) 

r collision time, sec 

r c collision time for multiple -reflection collision, sec 
t q collision time for single -reflection collision, sec 

transition time for ion-molecule charge transfer, sec 
cp azimuthal angle for translational motion, rad 
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APPENDIX B 


EQUATIONS OF MOTION AND VALUES OF MOLECULAR CONSTANTS 

The Lagrangian energy expression for the ion - linear molecule collision system 
depicted in figure 1 (a) is 


L = T - V = “ (X 2 + Y 2 + Z 2 ) + -1 (k 2 + p 2 sin 2 |) 




(Al) 


where m, I, p, a, and e are the reduced mass, moment of inertia, dipole moment, 
electronic polarizability, and electronic charge, respectively. The potential energy 
term is written in atomic physics units (i. e. , substitute in denominator to con- 

vert to SI units) . 

The equations of motion are simply 


d^ 

dt 



(A2) 


where q = X, Y, Z, and 77 . 

The translational coordinates equations are 


X = -M- (r 2 sin | cos 77 - 3XF’) - 2ae2x 
mr 5 mr 6 

Y = Me ~ (r 2 sin £ sin 77 - 3YF* ) - 2 . ae ^. Y 

5 v 6 

mr mr 


(A3) 


(A4) 


Z = 


~ y- e (r 2 cos £ - 3ZFp - 
mr 5 


2ae 2 Z 


mr 


where the potential function 


F^ = X sin £ cos 77 + Y sin £ sin 77 + Z cos £ 
24 


(A5) 



The corresponding rotational coordinate equations are 


\ = T^sin £ cos £ - /J,e (Z sin £ + X cos £ cos 77 - Y cos £ sin 77) 

v 3 


(A6) 


77 = 



(Y cos 77 - X sin 77) - 21^77^ cos | 


Ij sin £ 


(A 7 ) 


Numerical singularity difficulties caused by the sin £ term in the denominator of 
equation (A 7 ) prevent an accurate calculation of the azimuthal acceleration 77. This 
singularity is avoided by transforming to a second Euler angle system whose polar axis 
is orthogonal to that of the first system (refs. 3 , 8 , and 9 ). The form of the equations 
for symmetric top targets (see fig. 1 (b)) has been extensively discussed in reference 8 . 
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APPENDIX C 


INTEGRATION AND PLOTTING PROCEDURES 

The differential equations are integrated by using the variable-step-length Runge- 
Kutta scheme discussed in appendix C of reference 9. The error introduced in each in- 
tegration step is estimated by considering the difference between the theoretically con- 
served total energy T + V before and after the step. If this difference is greater than 
a specified tolerance (proportional to the step length), the step length is reduced and the 
integration procedure repeated. However, if the energy difference is less than the 
specified tolerance, the step is accepted and a new length is computed for the next 
step. 

The initial conditions were generated randomly in the manner described in refer- 
ences 3, 8, and 9. These initial conditions and values at the end of each successful in- 
tegration step are stored in a plotting array for each of the quantities to be displayed by 
computer plotting. After the integration is completed, the plotting subroutines are 
called. Each plotting array contains N + 1 values, where N is the number of success- 
ful integration steps. A typical range for N is approximately 500 to 1000. For each 
plot the contents of one array (e. g. , coordinate angles, relative velocity, or rotational 
energy) are taken as ordinate values to be plotted against the abscissa values (usually r, 
the ion -molecule separation). The consecutive points so determined are then connected 

by N straight lines. The total energy variation tolerated over the entire integration is 
4 -3 

usually set at 10 to 10 . Standard fourth -order Runge-Kutta and fifth and sixth orders 

are most often used. The resulting step lengths for these combinations are sufficiently 
short so as to give acceptably smooth curves. An exception to this regular behavior is 
observed for unusually low rotational energy cases. The trajectory solutions are also 
slightly sensitive to step -size histories (i. e. , the order in which the trajectories are 
solved). The number of multiple reflections may change with step-size history, but this 
is not felt to prejudice the statistical worth of the results. 

As noted in reference 9, this integration scheme has yielded more accurate results 

and significantly shorter computer run times than fixed-step-length routines. For ex- 

+ 4 

ample, the system Ar + CO with a relative translational velocity v of 5. 5x10 centi- 
meters per second and a rotational energy Ej^ « kT^ for T^ = 300 K required an 
average computer time per case of 10 to 15 seconds, including plotting. 

The plotting hardware and the subroutines required to operate it are described in 
reference 12. The DD280 recorder consists of a cathode- ray display tube and a 35- 
millimeter still camera. Calls to the plotting subroutine produce images on the display 
tube, which is an on-line output device to the Lewis IBM 360/67 computer. The 35- 
millimeter camera photographs the image and advances the film to the next frame. 
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After developing, the film can be studied in a viewer, printed on photographic paper, or 
submitted for reproduction in a publication. The plotting subroutines are quite flexible, 
permitting programmer-designated plot titles, axes labels, point or curve labels, grid 
selection, plot size and margin adjustment, plots extending over several frames, and 
multiple plots on a single frame. 
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